Epidemic dynamics in finite size scale-free networks 
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Many real networks present a bounded scale-free behavior with a connectivity cut-off due to phys- 
ical constraints or a finite network size. We study epidemic dynamics in bounded scale-free networks 
with soft and hard connectivity cut-offs. The finite size effects introduced by the cut-off induce an 
epidemic threshold that approaches zero at increasing sizes. The induced epidemic threshold is very 
small even at a relatively small cut-off, showing that the neglection of connectivity fluctuations in 
bounded scale-free networks leads to a strong over-estimation of the epidemic threshold. We pro- 
vide the expression for the infection prevalence and discuss its finite size corrections. The present 
work shows that the highly heterogeneous nature of scale-free networks does not allow the use of 
homogeneous approximations even for systems of a relatively small number of nodes. 
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In the last years it has been recognized that a large 
number of physical, biological, and social networks ex- 
hibits complex topological properties jl], g|. In particu- 
lar many real world networks show the small-world phe- 
nomenon, related to a very small average path length 
between nodes ||, |j| . More strikingly, in some cases this 
property is associated to a scale-free connectivity distri- 
bution, P(k) ~ fc~ 2-7 , with < 7 < 1, where k is the 
number of links connected to a node . This scale-free 
nature is associated to a large heterogeneity in the con- 
nectivity properties of the system. Since the second mo- 
ment of the connectivity distribution (fc 2 ) is diverging 
when increasing the network size, the connectivity fluc- 
tuations in scale-free (SF) networks do not have an in- 
trinsic bound and diverge in the infinite system size limit. 
Scale-free properties have been observed in many real sys- 
tems such as the Internet || ^, Q and the World-Wide- 
Web j|, ||, food- webs, protein, and neural networks ||. 
A very important example of scale-free networks is also 
found in the web of human sexual contacts ]To| . This is 
a particularly relevant case since the unambiguous defi- 
nition of contacts (links) is often missing in the analysis 
of social networks. 

Since the Internet and the web of human sexual con- 
tacts appear to be scale-free, the study of epidemics and 
disease dynamics on SF networks is a relevant theoreti- 
cal issue in the spreading of computer viruses and sexu- 
ally transmittable diseases. In heterogeneous networks, 
it is well-known that the epidemic threshold decreases 
with the standard deviation of the connectivity distri- 
bution jn], and this feature is amplified in SF net- 
works, which have diverging connectivity fluctuations in 
the limit of infinite network size. Indeed, it was first 
noted in Ref. that, in infinite SF networks, epidemic 
processes do not possess an epidemic threshold below 
which diseases cannot produce a major epidemic out- 
break or the inset of an endemic state. The absence of an 
intrinsic epidemic threshold has been found in both the 
susceptible-infected-susceptible (SIS) model | D| and the 
susceptible-infected- removed (SIR) model [|14 , |15| ] in infi- 



nite SF networks. The immunization policies are as well 
very much affected by the SF nature of the connectivity 
distribution @, 0. 

As customarily encountered in nonequilibrium statis- 
tical systems Q, it has also been pointed out that in 
finite systems an epidemic threshold is induced by finite 
size effects |Q. Real systems are actually made up by a 
finite number of individuals which is far from the ther- 
modynamic limit. This finite population introduces a 
maximum connectivity k c , depending on N , which has 
the effect of restoring a bound in the connectivity fluctu- 
ations, inducing in this way an effective nonzero thresh- 
old. More generally, we can consider a class of bounded 
scale-free (BSF) networks, in which the connectivity dis- 
tribution has the form P(k) ~ fc _2_7 /(/c/fc c ), where the 
function f(x) decreases very rapidly for x > 1 [^9|. The 
cut-off k c can be due to the finite size of the network or to 
the presence of constraints limiting the addition of new 
links in an otherwise infinite network M. 

In this paper we present an analytical study of the 
SIS model in BSF networks with a generic connectivity 
exponent 7 (0 < 7 < 1), focusing on the effects intro- 
duced by a finite cut-off k c . We analyze the case of a 
hard cut-off, f(x) = 6(1 — x), where 6(x) is the Heav- 
iside step function, as it happens in growing networks 
with a finite number of elements. We consider as well a 
soft exponential cut-off, f(x) = exp(— x), as often found 
in systems where physical constraints are at play. We 
derive the behavior of the epidemic threshold as a func- 
tion of k c and the network size N, and find that even for 
relatively small networks the induced epidemic thresh- 
old is much smaller than the epidemic threshold found in 
homogeneous systems. This confirms that the SF nature 
cannot be neglected in the practical estimates of epidemic 
and immunization thresholds in real networks. We also 
provide the explicit analytic form for the epidemic preva- 
lence (density of infected individuals) in BSF networks. 
The results presented here can be readily extended to the 
SIR case. 

In order to estimate the effect of k c in epidemics 
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on BSF networks we will investigate the standard SIS 
model 1 20 . This model relies on a coarse-grained descrip- 



tion of individuals in the population. Namely, each node 
of the graph represents an individual and each link is a 
connection along which the infection can spread. Each 
susceptible (healthy) node is infected with rate v if it is 
connected to one or more infected nodes. Infected nodes 
are cured and become again susceptible with rate <5, defin- 
ing an effective spreading rate A = v/8 (without lack of 
generality, we set 8 = 1). The SIS model does not take 
into account the possibility of individuals removal due to 
death or acquired immunization p(| , and individuals run 
stochastically through the cycle susceptible — » infected — ► 
susceptible. This model is generally used to study infec- 
tions leading to endemic states with a stationary average 
density of infected individuals. In order to take into ac- 
count the heterogeneity of SF networks, we have to relax 
the homogeneity assumption used in regular networks, 
and consider the relative density pk(t) of infected nodes 
with given connectivity fc; i.e., the probability that a node 
with k links is infected [pL3|| . The dynamical mean- field 
equations can thus be written as 



dpk(t) 
dt 



-p k (t) + Xk [l-p k (t)]Q(p(t)). 



(1) 



The first term in Eq. (|IJ) considers infected nodes be- 
coming healthy with unit rate. The second term rep- 
resents the average density of newly generated infected 
nodes that is proportional to the infection spreading rate 
A and the probability that a node with k links is healthy 
[1 — Pk(t)] and gets the infection via a connected node. 
The rate of this last event is given by the probability 
<d(p(t)) that any given link points to an infected node, 
which has the expression |13| 



e(p(t)) = (k)- 1 J2kP(k)Pk(t). 



(2) 



By solving Eqs. (0) and (j|) in the stationary state 
[dpk{t)]dt = 0] we obtain the self-consistency equa- 
tion 



(3) 



where is now a function of A alone. The self-consistency 
Eq. (^) allows a solution with 0^0 and pk ^ only if 
the condition A (fc 2 ) / (fc) > 1 is fulfilled fig] , defining the 
epidemic threshold 



A c = 



(fc 2 )' 



(4) 



In other words, if the value of A is above the threshold, 
A > A c , the infection spreads and becomes endemic. Be- 
low it, A < A c , the infection dies out exponentially fast. 
This result implies that in infinite SF networks with con- 
nectivity exponent < 7 < 1, for which (fc 2 ) — > 00, 
we have A c = 0. This fact implies in turn that for 



any positive value of A the infection can pervade the 
system with a finite prevalence, in a sufficiently large 
network [Fl3|. While this results is valid for infinite SF 
networks,^!; 2 ) assumes a finite value in BSF networks, 
defining an effective nonzero threshold due to finite size 
effects as usually encountered in nonequilibrium phase 
transitions This epidemic threshold, however, is 

not an intrinsic quantity as in homogeneous systems and 
it vanishes for a increasing network size or connectivity 
cut-off. In order to calculate the precise effects of a fi- 
nite k c we consider two different cases of connectivity 
cut-off. At first instance we consider a soft exponential 
cut-off with characteristic connectivity k c . This case cor- 
responds to those real networks in which external factors 
set up an upper limit to the connectivity 0. The net- 
work can have an infinite number of elements but the 
power-law connectivity distribution decays exponentially 
for large values of A:. In order to perform explicit calcula- 
tions we use a continuous approximation that substitutes 
the connectivity by a real variable k in the range [to, 00), 
where m is the minimum connectivity of the network. 
The connectivity probability distribution in this case is 
P(k) = Ak^ 2 ^ 1 exp(— k/k c ), where A is a normalization 
factor. The effective nonzero epidemic threshold A c (fc c ) 
induced by the exponential cut-off is given by 

£k-^cM-k/k c )dk 
( )_ C k-^M-k/k c )dk ' (5) 

which, after integration, yields 



A c (fc c ) = k 



-1 r(-7,m/fc e ) 
r(l-7, m/k c y 



(6) 



where T(x, y) is the incomplete Gamma function fl21| . For 
large k c we can perform a Taylor expansion and retain 
only the leading term, obtaining for any < 7 < 1 



A c (fc c 



1 



m7r(l — 7) 



(k c /m 



,7-1 



(7) 



The limit 7 — » 1 in Eq. (||) corresponds to a loga- 
rithmic divergence, yielding at leading order A c (fc c ) ~ 
(to ln(A: c /m)) _1 . In all cases we have that the epidemic 
threshold vanishes when increasing the characteristic cut- 
off. For large fc c , the average connectivity is virtually 
fixed and given by (k) = (7 + 1)771/7, f° r an y 7 > 0. 
It is interesting thus to compare the intrinsic epidemic 
threshold obtained in homogeneous networks with negli- 
gible fluctuations and the nonzero effective threshold of 
BSF networks. The intrinsic epidemic threshold of ho- 
mogeneous networks with constant node connectivity (k) 
is given by A" = (fc) -1 @ If we compare BSF and 
homogeneous networks with the same average connectiv- 
ity (fc) = (7 + 1)771/7 we obtain that the ratio between 
the epidemic thresholds is given by 



A c (fc c ) (7 + I) 

ah - 7 2r(i- 7 ) 



(kc/m) 



7-1 



(8) 
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FIG. 1: Ratio between the effective epidemic threshold in BSF 
networks with a soft exponentiai cut-off k c and the intrinsic 
epidemic threshold of homogeneous networks with the same 
average connectivity, for different values of 7. 



This clearly shows that even in the case of a connectivity 
cut-off the effective epidemic threshold in BSF networks 
is much smaller than the intrinsic threshold obtained in 
regular networks. In Fig. [l] we plot the ratio obtained by 
using the full expression for X c (k c ), Eq. (^). It is striking 
to observe that, even with relatively small cut-offs (k c ~ 
10 2 — 10 3 ), for 7 ss 0.5 the effective epidemic threshold of 
BSF networks is smaller by a factor close to 1/10 than the 
intrinsic threshold obtained on homogeneous networks. 

As a second kind of finite size effect we consider the 
presence of a hard cut-off k c . Since SF networks are of- 
ten dynamically growing networks, this case represent a 
network which has grown up to a finite number of nodes 
N. The maximum connectivity k c of any node is related 
to the network age, measured as the number of nodes N, 
by the scaling relation pi 



(9) 



where m is the minimum connectivity of the network. 

In this case the network does not possess any node 
with connectivity k larger than fc c , and we can think in 
terms of a hard cut-off. Using again the continuous k 
approximation, the normalized connectivity distribution 
has now the form 



P(k) 



(1 + 7)771 



1+7 



1 - (k c /m) 



-1-7 



k- 2 -^6{k c - k), 



(10) 



where 9{x) is the Heaviside step function. The finite 
size induced epidemic threshold A c (fc c ) is given by the 
expression 



\ c (k c 



r° k-^dk 

J m 

f kc k~^dk 



(ii) 



Evaluating the above expression we obtain at leading or- 
der in k c /m: 



A c (/c c ) 



7 



7m 



(k c /m 



(12) 



In this case the hard cut-off k c can be expressed as a 
function of the network size N by using the scaling re- 
lation Eq. (||) and we can obtain the effective epidemic 
threshold as 



\c(N) 



1 



7771 



(13) 



This expression is valid for any < 7 < 1 , while for 7 = 1 
we obtain at the leading order the logarithmic behavior 
\ C (N) ~ 2(to ln(iV)) _1 . Also in this case we have that 
the effective epidemic threshold is approaching zero for 
increasing network sizes, and it is worth comparing its 
magnitude with the corresponding intrinsic threshold in 
homogeneous networks with identical average connectiv- 
ity. In Fig. U we report the ratio A c ( N) / A^ for different 
sizes of the SF network. It is striking to notice that for 
7 = 0.5, small networks with N ~ 10 4 exhibit a finite 
size induced epidemic threshold that is close to be one 
order of magnitude smaller than the intrinsic epidemic 
threshold of a homogeneous network. 

In order to find the prevalence behavior we have to 
solve Eq. (|^) in the continuous approximation, 



e = (fc)- 1 



\ ( \ V> 

tp ^TTxek dk - < 14) 



and use the value of <d to compute the density of infected 
sites p as 



P = Y,PkP[k)= / P(k)- 



xek 



xek 



-dk, 



(15) 



where P(k) is given by Eq. (|l0|). In the absence of any 
cut-off (k c — ► 00) and in the thermodynamic limit (JV — » 
00) the prevalence scales as p ~ A 1 ^ 1-7 ' if < 7 < 1, 
and as p ~ exp(— I/777A) if 7 — 1 Jl3[ . Accordingly with 
the absence of the epidemic threshold, the prevalence is 
null only if the spreading rate is A = 0. In the case of a 
hard cut-off we can integrate Eq. (|l4|), neglecting terms 
of order (k c /m)~~" in the P(k) distribution, to obtain: 



e 



7m 7 AG 



fc" 7 



-dk 



1 + A6fc 
= F(l,7,l + 7 ,-[Ae777]- 1 ) 

-(fc c /m)- 7 F(l, 7 , 1 + 7, -[A9A:,]- 1 ), 

where F is the Gauss hyper-geometric function plj . For 
a fixed k c one can expand both hyper-geometric functions 
on the right hand side in the previous equation, keeping 
the most relevant terms in O and considering afterwards 
the limit of large k c . The final solution for O is then 
given, at leading order in (fc c /m), by 



e 



1 2 



7 



mA 1 — 7 



A - 



1-7 
7777 



(16) 



By evaluating the integral in Eq. (|T^) and keeping the 
leading term in and k c we finally obtain the infection 
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FIG. 2: Ratio between the effective epidemic threshoid in 
BSF networks with finite size N and the intrinsic epidemic 
threshoid of homogeneous networks with the same average 
connectivity, for differen values of 7. 



prevalence as 

_ ( 7 +l)(2- 7 ) (k L 



A7(l — 7) \m 



A - 



1-7 (h 

7m 



7-1 



Inserting the scaling relation Eq. (^|) between the maxi- 
mum connectivity k c and the network size N we are led 
to the final expression 



^-1/(7+1) (A - X C (N)). 



(17) 



That is, the finite size of the network induces a standard 
mean-field transition at the induced epidemic thresh- 
old A c (AT), given by Eq. (|l3|). As can be seen from 



Eq. (KL7p, however, the prevalence is depressed by a factor 
N- 1 /Vi+ 1 ) f rom the corresponding value for a homoge- 
neous network. The above calculations can be repeated 
along similar lines in the case of a soft exponential cut-off, 
obtaining similar results. 

It is worth remarking that similar results hold as well 
for the SIR model. Despite this model confers perma- 
nent immunity and does not allow for a stationary state, 
the epidemic threshold over which an epidemic outbreak 
occurs has the same analytic form A c = (k) / (k 2 ^ Jl5| . 
Thus, the present results for the effect of finite size and 
the induced epidemic threshold can be readily exported 
to the SIR case. The calculation of the epidemic preva- 
lence is different due to the different evolution equations, 
but recovers the same onset of an induced mean-field 
transition at the effective threshold \ C (N). Interest- 
ingly, similar results have been obtained in the analysis of 
the resilience to damage of finite size scale-free networks 

nn 

In conclusion, we have shown that the SF networks 
weakness to epidemic agents is also present in finite size 
networks. Using the homogeneity assumption in the case 
of SF networks will lead to a serious over-estimate of the 
epidemic threshold even for relatively small networks. 
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